IEEE Xplore Citation 



http://iecexplore.ieee.org/search97/s97i... f=&Col I Std=&sortfield=pyr&sortorder=desc 



Help FAQ Terms Release Notes Rele 



Welcome to SE££Ai?i0r 



Q" Home 
Q- Log-out 



Q- Journals 
& Magazines 

0- Conference 
Proceedings 

0~ Standards 



O" By Author 
0~ Basic 
O" Advanced 

O-Joiftl&BE 

O Establish 3EEE 
Wet) Accost 

r3 Print Formal 



SEARCH RESULTS fPDF Full-Text (260 KB^I PREVIOUS NEXT 

New Fast GPS code-acquisition technique using FFT 

- Van Nee, D.J.R.; Coenen, AJ.R.M. 

Delft Univ. of Technol., Netherlands 

This paper appears in: Electronics Letters 

On page(s): 158 - 160 

17 Jan. 1991 

Volume: 27 Issue: 2 

ISSN: 0013-5194 

References Cited: 4 

CODEN: ELLEAK 

INSPEC Accession Number: 3846752 



Abstract: 

A new spread-spectrum code-acquisition technique for the navigation systems 
Navstar/GPS and Glonass is introduced. This technique uses the FFT to compu 
correlation function, thereby eliminating the time-consuming code phase shift p 
Comparisons with existing systems show a theoretical reduction in acquisition t 
about 2000 times. 
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length ratios of u/L, v/L and w/L are 0-3, 0*4, and 0 3, respec- 
tively. As indicated in Fig. 2, more than -20dB crosstalk 
level of bar slate can be expected over the range from 3 to 4 of 
t$ T Jb0 Tl . By designing the length of each section (u, i\ w) 
adequately, a voltage controlled optical power splitter can be 
produced as shown in Fig, 3. Uniform &0 elements, whose 
Lfl ry and L// re were about I and 3, were used to measure the 
ratio a. The measured ratio a was from 3 6 to 3-8 in our 
experiment and was bigger than the previously reported 
values. 3,5 This is presumably due to the difference in mode 
confinement conditions for TM and TE modes between mea- 
sured directional couplers. However, a of around 3-7 is 
favourable for the performance of the device proposed. 




Fig. 3 Calculated characteristics as optical power splitter 
u/L - 0 3675; o/L - 0 265; w/L - 03675; a = 3-6 

Fabrication: The coupling length conditions for the TM and 
TE modes were satisfied when 750A-thick titanium, 7^m 
width, and 7 fim spacing patterns were diffused at 1050'C for 
6 h in a wet Ar atmosphere. The length of each section was 
determined based on the measured Afi TM /\f3 Tt . The coupling 
effects of the curved guide regions (A - 40 mm) at the both 
sides of directionaJ couplers for the TM and TE modes were 
also considered to define electrode photomask patterns. These 
coupling effects were measured by comparing the output light 
intensities from the guides of directional couplers with theo- 
retical values, and their estimated values in terms of coupling 
length were less than I mm for both modes. From these funda- 
mental data, devices with the length ratios of u/L « 0-375, 
v/L - 0-25, and w/L = 0-375 (L = 22 mm) were fabricated for 
optical power splitters. 

Experimental results: The device characteristics were mea- 
sured by PMF endfire coupling at 1-31 fim wavelength. A 
computer controlled measurement system with a IR camera 
and a voltage generator was used. Digital video signals from 




applied voltage, V 



Fig. 4 Measured and calculated characteristics 

O measured for TM mode 
Q measured tor TE mode 

TM mode 

TE mode 

u/L = 0-375; c/L = 0-25; w/L = 0-375; a = 3-5; L= 22 mm 



the IR camera were integrated and computed to measure 
crosstalk values. 

Fig. 4 shows the measured optical output intensities for 
both polarisations as functions of the applied voltages. In the 
Figure, calculations are also plotted. The ratios of the coup- 
ling length to the complete coupling lengths were about 0-9 
for the TM mode and about 2-9 for the TE mode. As predict- 
ed by the calculations, the switching voltages for the TM and 
TE modes are coincident, and the applied voltage dependence 
of the output light intensities for the T M mode is also coin- 
cident with that for the TE mode. The switching voltage is 
relatively low, about 15 V. 

Conclusion: New polarisation-insensitive devices, based on the 
modified configuration of three section alternating have 
been produced in Z-cut LiNbOj. Polarisation-insensitive 
switching and power splitting characteristics have been 
obtained 

N. KWZUTA 26th November 1990 

K. TAK.AKURA 

Skunadiu Co. Opto-electronia Department 
23B5-I) Ihiyama Atsuai, Kanagawa 243-02, Japan 
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NEW FAST GPS CODE-ACQUISITION 
TECHNIQUE USING FFT 



Indexing terms: Correlation, Fast Fourier transforms 

A new spread-spectrum code -acquisition technique Tor the 
navigation systems Navstar/GPS and GlotuusU introduced. 
This technique uses (he FFT to compute the corrclaiiun 
function, thereby eliminating ihe (ime-coosuming code phase 
shift process. Comparisons with existing systems show a 
theoretical reduction in acquisition lime of about 2000 limes. 



Introduction: A short acquisition lime is very important for a 
standard positioning service Navstar/GPS or Glonass recei- 
ver, especially in an urban environment where satellites are 
often visible (or a few seconds on!) In this case, an acquisition 
time of less than a few seconds is desired to avoid the receiver 
having to continuously remain in the acquisition phase. We 
discuss the code acquisition time only, that is. the time needed 
to align the incoming code and the local code within one chip. 
The proposed new acquisition technique is a result of our 
research into satellite navigation in an urban environment. 

Noncoherent correlator: The most frequently used code acqui- 
sition system is the noncoherent correlator, shown in Fig. l.'- : 
The incoming signal x{i) consists ol noise plus the GPS signals 
which have a carrier frequency of I 6 GHz, and are binary 
phase modulated by a Gold sequence of 1023 chips with a 
chip rate of I MHz, and by a 50 Hz daia stream. This signal is 
converted to baseband and coherently correlated with the 
local code for NLT t seconds. Here L is the code length 
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(= 1023 chips for the C/A code). T t is the chip time aod N is 
an integer £l. This time is chosen short enough to ensure 
thai the presence of data and carrier Doppler shift will noi 
cause a great degradation in performance. 3 Next, K sequential 
correlations are noncoherently summed 10 produce one corre- 
lation point at a sufficiently high signal-to- noise ratio. The 
total integration time KNLT t is called the dwell time. 



ReU(nl) 



sm(wt) 



lmU(n)l 



Flj. 1 Noncoherent correlator in tine domain 

To achieve acquisition with the noncoherent correlator, the 
above procedure has to be performed for all possible code 
phases and also for a number of possible Doppler frequencies, 
because the carrier phase has to remain nearly constant 
during individual integration times NLT t . As a result, the 
acquisition time is proportional to the number of cells, which 
is defined as the product of the number of code phase steps 
and the number of carrier frequency steps. For GPS. a typical 
value of the number of cells is 20460 for haJf a chip phase step 
(2046 steps for the C/A code) and I k Hz carrier frequency step 
at a correlation time of 1 ms {N - 1) and a carrier Doppler 
range of 10 kHz, neglecting the contribution due to the speed 
of the receiver. 

Parallel search techniques using the FFT: To reduce the acqui- 
sition time, cells can be searched in parallel by taking the FFT 
of the complex samples at the points I and 0 in Fig. I. 3 At the 
moment that the incoming code and the local code have the 
same phase, a carrier component will be present in the spec- 
trum, which is visualised by the FFT. Using this technique 
only 2046 phase steps remain, reducing the acquisition time 
ten times in comparison with the first system, assuming the 
FFT can be computed within the dwell time. 

In the above system, a parallel search in frequency is per- 
formed to eliminate 10 frequency search steps. In our pro- 
posed system, a parallel search in time is performed, i.e. all 
points of the correlation function are calculated from the 
same input sequence, with a theoretical gain in acquisition 
time of 2046 times in comparison with the conventional cor- 
relator, which is much more than the gain of the previous 
system that also uses a digital signal processor, but in a com- 
pletely different way. 
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Fit. 2 Noncoherent correlator via frequency domain 

The correlation operations are implemented as follows (see 
Fig. 2). A digital signal processor loads 2NL complex samples 
with a lime spacing of half a chip and then performs the 
correlation operation 



sl ■■ t 



form* 0,1 2L-I (1) 



The processor stores all 2NI. points of R J (m) in an accumula- 
tor. This process is repeated until K separate correlation func- 
tions have been summed. 

Unfortunately, a straightforward calculation of R(m) 
requires a large number of operations, which is proportional 
to {NL) 2 . Much computing time can be saved, however, if the 
correlation function is calculated via the frequency domain: 4 



where • represents convolution. X(k) is the spectrum of xfn), 
f\k) is the complex conjugate of the spectrum of p{n + m) and 
F~ 1 is the inverse Fourier transform. 

The fastest FFT algorithms require the number of points to 
be radix 2. We therefore add two zeros in both incoming and 
local code with a spacing of half the code length. In fact, we 
slightly deform the codes, which yields the advantage of fast 
processing but has the disadvantage of a somewhat larger 
crosscorrelation. The latter effect, however, is negligible com- 
pared to the thermal noise level in the acquisition phase. 

Because the digital signal processor computes the entire 
correlation function in one dwell lime, this technique is theo- 
retically 2046 (2L) times faster than the previously mentioned 
parallel frequency search method, which uses the same 
amount of hardware. In practice, this gain will only be 
achieved if the signal processor is fast enough to compute 2K 
times a 2NL point FFT within one dwell time, Also, this has 
to be done for at least ten difTcreni carrier frequency steps; for 
each frequency step the input signal is frequency shifted 
through a software-implemented inuge-rejection mixer. For 
the typical values K * 20. N = 1 and ten frequency steps, this 
would mean that ten 2048 point FFTs have to be computed 
within 500 /is. Such signal processors actually exist but are still 
extremely expensive. However, even with the TMS320C30 
processor we use, an estimated minimum acquisition time of 
less than 500 ms per frequency step is possible, which then 
yields a gain in acquisition time of eight times compared to 
the parallel frequency search system. In the case of reacquisi- 
tion, which will occur frequently in an urban environment, the 
carrier frequency will often be known accurately enough to 
reduce the number of frequency search steps to only one, so 
the gain in acquisition time of the TMS320C30 can be 
increased by a factor ten. 




Fig. 3 Normalised mean acquisition lime as junction of the sianal-to- 
noise ratio for K = 20 and K » 4 

K - 20 

K -4 

+ simulated results for K - 20 
x simulated results for K » 4 
L = 1023 c * 5 

Finally, to demonstrate that the system actually works, Fig. 
3 shows a plot of the mean acquisition time as a function of 
the signal-to-noise ratio E t )N n which is calculated using the 
equations described in Reference 3. For a number of difTcreni 
signal-to-noise ratios, we simulated a OPS spread-spectrum 
signal with additive Gaussian noise. We performed 100 simu- 
lations for each value of the signaJ-to-noisc ratio, and counted 
the number of successful acquisitions, which is a measure of 
the acquisition probability p that can be substituted in the 
equation for the mean acquisition time in the case of a single 
dwell acquisition procedure (for 
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The constant c . KLT t is the penalty time in the case of a false 
alarm. As can be seen, the simulated plots have a close resem- 
blance to the calculated curves. The shift of about half a dB 
for K = 20 is most probably due to the approximations in the 
calculation, so wc can conclude that the performance of the 
systems in Fig. I and 2 are the same, except for a constant 
factor in the acquisition lime. It should be noticed that the 
mean acquisition time in Fig. 3 ts multiplied by a certain 
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factor if the signal processor used is not processing in real 
time. 

Conclusions: By using a digital signal processor in a hitherto 
unusual way. the code acquisition time or a GPS receiver can 
be reduced considerably. The fastest performance is achieved 
when the incoming code and the local code are slightly 
deformed by adding zeros, allowing the use of the FFT. 

Acknowledgment: The authors would like to thank Mr. E. 
Theunissen and Mr. T. Lamerigts for their work on the 
TMS320C30 processor and Prof. D, van Willigen and Prof. R. 
Prasad for their helpful comments. 
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A. J. R. M. COENEN 
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INJECTION- LOCKING OF G-SWITCHED 
AIGaAs LASER WITH FAST SATURABLE 
ABSORBER 

Indexing terms: Semiconductor lasers. Lasers 

Injection of a weak CW lam beam in a Q-switcbed AIGaAs 
laser diode with a fesi saturable absorber is shown to 
produce powerful single-mode picosecond pulses at 0-82 >mv 
The saturable absorber regions are obtained by deep implan- 
tation of heavy ions through the diode facets. Single-mode 
operation is achieved with CW injection powers as low as 
50jiW. Peak powers exceeding I S W are detected at the laser 
output. A time- resolved spectroscopy of the laser pulses 
reveals an overall do vnchirp of l-Snm. 



Introduction: Q-s witching of (Al)GaAs laser diodes is now rec- 
ognised as being a very efficient and simple technique to 
produce picosecond purses in the 0-8-0- 85 /an region. 1 "* 
Recently, 5ps pulses with peak powers up to 5W have been 
observed by using fast saturable absorbers. 4 3 Pulses produced 
by this technique, however, generally exhibit poor spectral 
characteristics that are not desirable for a number of applica- 
tions. We demonstrate the possibility of controlling the laser 
spectrum with the injection of a weak CW beam into the 
cavity. Owing to the use or a very fast saturable absorber, 
output powers are ten limes higher than those reported in the 
previous works on injection -locking of pulsed semiconductor 
lasers. 6 Moreover, lime-resolved spectroscopy of the single- 
mode laser pulses reveals a downchirp of very large ampli- 
tude, which enables us to predict the feasibility of 
bandwidth-limited subpicosccond pulses. 

Experimental setup: The experimental setup is shown in Fig. 1. 
The semiconductor laser is a Ga As/ AIGaAs gain-guided 
doublc-hctcrostructurc already described in previous works. 1 3 
The active region is 200 long. The emission wavelength is 
in the 820-830 nm range. Internal 0 -switching is obtained by 
producing regions of a saturable absorber near the mirrors of 
the laser resonator by ISMeV implantation of nitrogen ions. 
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With electrical pump pulses of 40 V amplitude, 1-5 ns riselime 
and 100 kHz repetition rate, the laser emits multimode optical 
pulses of 5-1 Ops full width at half maximum (FWHM). The 
energy per facet is 1 5 pJ. 



f-V r freer (ter 




fi sapphire 1 
laser 



Fif. 1 Experimental setup 

The laser mode control is presently achieved by injecting a 
weak CW tunable Ti : sapphire laser beam through one facet 
of the laser diode. A 0-1-NA /= 20 mm lens is used to focus 
the injection beam in the resonator. The CW infrared power 
incident onto the sample typically ranges from 10 to 100>W, 
only 10% of which is estimated to be injected into the trans- 
verse resonator mode. The residual amount of power trans- 
mitted from the opposite facet of the laser diode is still one 
order of magnitude lower and typically represents less than 
10% of the average output power delivered by the diode. 

For spectral measurements, the output beam is diffracted 
with a I200g/mm grating after being collimated with a 
0-35-NA J = 1 0mm microscope objective. An /= 50 cm lens 
is then used to focus the diffracted beam either onto an infra- 
red camera or a fast Si avalanche photodiodc. The camera is 
used for rapid control of the laser spectrum. Injection- locking 
is seen to occur when one of the laser spots becomes much 
brighter than the others. For a more precise analysis, the 
camera is replaced by the fast photodiode, the latter having a 
detecting aperture of 50 /un. A spectral resolution of about 1 A 
is obtained when scanning the photodetector across the beam. 

The same arrangement is used to qualitatively analyse the 
wavelength evolution during the pulse. As shown in Fig. 1, the 
photodetector signal is fed into the 7S12 sampling unit of a 
Tektronix 7854 oscilloscope equipped with the S4 sampling 
head. Using the manual mode of operation, the signal is 
analysed in a temporal window of — 2Ups, the position of 
which can be adjusted with a precision of a few picoseconds 
with respect to the triggering pulse. The 7S12 unit output is 
then recorded while slowly scanning the phutodeteclor across 
the diffracted beam. The operation is repeated for different 
positions of the temporal window. Because of the longer 
lOOps response time of the photodetector as compared to the 
^10ps width (FWHM) of the optical pulses, the photo- 
detector actually acts as an integrator. In other words, the 
procedure described above allows us to measure the time- 
averaged spectrum of the laser pulse for various averaging 
times in the tens of picoseconds scale. Although the resolution 
is limited by the width of the temporal window, we are able to 
determine the sign and the total amplitude of the wavelength 
shift during the pulse. 

Experimental results: The efficiency or the injection-locking 
process is illustrated in Fig. 2. The bottom curve represents 
the multimode spectrum of the laser without injection. The 
pulse energy spreads over a dozen of the longitudinal modes, 
which approximately corresponds to a spectral width of S nm. 
The asymmetrical shape of the spectrum envelope indicates 
the presence of chirping effects. The top curve in Fig. 2 rep- 
resents the laser spectrum under injection-locking conditions. 
Most of the pulse energy is now concentrated within a narrow 
spectral region located on the long wavelength side of the 
injection. The width of the spectral region essentially depends 
on the electrical pumping and typically ranges from 12 to 
l -6nm. The fraction of pulse energy contained within this 
region depends on the injection parameters. It increases with 
the injection power. At a given injection power, an optimum 
value is found when the injection is tuned near the gain 
maximum of the laseT diode. For the case of Fig. 2, the injec- 
tion power is 50jjW, the injection laser is tuned at 822 nm and 
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